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A.  Overview 

A  key  aspect  of  the  captioned  grant  has  been  the  simulation  of  photoelectron  spectra  of 
molecules  that  are  relevant  to  the  high  energy  density  materials  program  and  involve  states 
strongly  coupled  by  conical  intersections.  Our  AFOSR  funded  work  has  included  both  the 
development  of  new  computational  tools  to  determine  photoelectron  spectra  within  the  time 
independent  multimode  vibronic  coupling  model,  and  the  use  of  those  tools,  and  techniques 
developed  as  part  of  other  funded  research,  to  determine  photoelectron  spectra  for  the  selected 
azolides,  anion  precursors  of  the  heterocycle  radicals,  (CH)mN5  m,  for  m  =  0-4.  These  anions 
have  been  suggested  as  components  of  ionic  liquids  and  have  been  studied  experimentally  and 
computationally  in  Lineberger's  group.121,4  Key  to  our  ability  to  determine  photoelectron 
spectra  is  an  open-ended,  fine-grained  parallel  algorithm  for  solving  the  nuclear  Schrodinger 
equation  using  large  vibronic  bases.5  This  algorithm  (see  Sec.  B  below)  was  used  to  determine 
the  negative  ion  photodetachment  spectrum  of  the  anion  of  pyrazolyl,6  (CH)mN5.m,  for  m  =3, 
whose  determination,  as  recently  as  2006,  was  declared  computationally  unfeasible!4  Work  is 
currently  finishing  up  on  pyrrolyl,  tn  =  4,  is  significantly  in  progress  on  triazolyl  m  =  2  and  has 
begun  on  pentazolyl,  m  ~  0.  We  have  also  studied  a  novel  conical  intersection  topography  in 
excited  states  of  N3+  a  potential  photoionization  probe  of  cyclic  N3,  in  which  3  symmetry 
equivalent  seams  of  conical  intersection  are  found  in  close  proximity  to  a  symmetry  required 
seam.  Finally,  we  have  and  will  continue  to  do  work  on  the  role  of  conical  intersections  in  the 
photodissociation  of  nitramines,  which  has  been  suggested  as  being  relevant  to  the  ignition 
process  in  HMX,  RDX  and  related  species. 8,9 

B.  Computational  Tools 

(i)  Parallel  Lanczos  Algorithm 

M.  S.  Schuurman,  R.  A.  Young,  and  D.  R.  Yarkony,  Chem.  Phys.  347,  57  (2008). 
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We  developed  a  fine-grained  parallel  algorithm  to  solve  the  secular  equation  that  arises  in 
the  time  independent  multimode  expansion  of  the  vibronic  coupling  problem.50  The 
implementation,  which  defines  the  state  of  the  art  in  this  area,  can  handle  expansions  of  arbitrary 
length,  with  the  open-ended  character  of  the  algorithm  achieved  through  the  use  of  fine  grained 
parallelism  to  partition  the  trial  vectors.  The  algorithm  has  been  used  to  treat  vibronic 
expansions  comprised  of  over  1  billion  basis  functions.5 

C.  Photoelectron  Spectroscopy  of  the  Azolyls  (CH)mN5.m:  States  Strongly  Coupled  by 
Conical  Intersection 

These  studies  were  motivated  by  the  photoelectron  spectra  measured  in  Carl  Lineberger's 
laboratory,  which  we  initially  learned  about  at  an  AFOSR  contractor’s  meeting  in  2003.  We 
originally  showed  that  these  radicals  exhibited  low-lying  seams  of  conical  intersections  of  two 
and  three  electronic  states."  The  goal  of  the  captioned  research  effort  was  to  understand  what 
the  effect  of  this  complicated  series  of  conical  intersection  seams  has  on  the  photoelectron 
spectra.  The  primary  challenges  arise  from  the  strong  coupling  of  (at  least)  three  low-lying 
electronic  states,  and  the  consequent  low-lying  seams  of  two  and  three  state  conical  intersections, 
as  well  as  the  significant  number  of  internal  degrees  of  freedom  that  must  be  treated.  Our  recent 
methodological  advances  have  made  the  study  of  this  difficult  problem  tractable  using  the  time- 
independent  multimode  vibronic  coupling  method. 

(i)  The  photoelectron  spectrum  of  Pyrazolyl 

M.  S.  Schuurman  and  D.  R.  Yarkony,  J.  Chem.  Phys  127,064304  (2008). 

The  spectrum  of  the  pyrazolyl  radical,  (CH)3N2,  obtained  from  negative  ion 
photodetachment  of  an  electron  from  the  pyrazolide  anion,  proved  to  be  a  particularly  demanding 
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computational  problem.  Indeed,  as  recently  as  2006  the  determination  of  this  spectrum  was 
declared  computationally  unfeasible!4 

The  first  step  in  the  determination  of  the  photoelectron  spectrum  was  the  development  a 
quasi-diabatic  Hamiltonian,  Hd,  which  accurately  describes  the  electronic  structure  aspects  of  the 
problem  as  determined  by  ab  initio  multireference  configuration  interaction  wave  functions.  Not 
only  must  the  equilibrium  geometries  and  harmonic  frequencies  of  the  nominal  X2A2  ground 

state  and  low-lying  A  B ,  excited  state,  be  well  described  but  so  must  the  minimum  energy 

points  on  the  lowest  two-state  ( X2A2 ,  A  B{)  and  three-state  ( X2 A2 ,  A2 Bx ,  B2B2)  seams  of 

conical  intersection.  Hd  is  then  used  to  form  a  vibronic  Hamiltonian  in  the  multimode  vibronic 
coupling  approximation  in  a  basis  designed  to  limit  the  size  of  the  vibronic  expansion.  The 
vibronic  Hamiltonian,  of  dimension  greater  than  500  million,  was  diagonalized  to  obtain  the 
negative-ion  photoelectron  spectra  for  pyrazolide-/zv  Our  predicted  spectrum,  compares  quite 
favorably  with  the  measured  photoelectron  spectrum.4  See  Fig.  1  below. 
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Figure  1  :Expenmental(dots)  and  calculated  (solid  lines)  photoelectron  spectrum  of  pyrazolide-/o 
(ii)  The  photoelectron  spectrum  of  Triazolyl 
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The  photoelectron  spectrum  of  the  anion  of  triazolyl,  triazolide,  is  quite  challenging  to 
determine  computationally.  In  triazolyl,  the  combination  of  3  nitrogen  lone  pairs  and  a  half 
filled  5  electron,  5  orbital  orbital  conjugated  7t-system,  requires  that  four  low-lying  electronic 
states  be  considered.  The  proper  description  of  these  four  states  and  their  interactions,  which  is 
essential  for  the  determination  of  the  photoelectron  spectrum,  is  made  difficult  by  the  significant 
differences  in  the  equilibrium  geometries  of  the  four  states.  To  accommodate  this  situation  we 
are  employing  an  extension  of  the  procedure  used  to  construct  the  Hd  for  pyrazolyl,  a  procedure 
based  on  a  unique  normal  equations  approach.1  The  normal  equations  approach  allows  us  to 
incorporate  data  averaging  over  a  range  of  nuclear  geometries  while  assuring  quasi-diabaticity. 

Our  treatment  of  the  electronic  structure  of  triazolyl  is  quite  ambitious  being  an  1  1 
electrons  in  8  orbitals  active  space.  The  second  order  configuration  interaction  wave  functions, 
with  only  5  core  orbitals  kept  doubly  occupied,  comprised  more  than  200  million  configuration 
state  functions. 

At  present  we  are  in  the  process  of  determining  an  Hd  which  is  converged  with  respect  to 
the  normal  equations  averaging  procedure.  In  our  normal  equations  approach  for  triazolyl,  Hd  is 
constructed  using  data  from  both  the  region  of  the  ground  state  minimum  and  that  of  minimum 
energy  crossing  point.  Figure  2  below  compares  the  experimental  photoelectron  spectrum 
(upper  panel)  to  our  simulations  (lower  panels),  which  currently  employ  over  200  million 
vibronic  basis  functions.  The  two  lower  panels  juxtapose  spectra  obtained  with  the  minimum 
energy  crossing  (  lower  left)  and  ground  state  energy  minimum  (lower  right)  as  the  origin  for 
Hd,  respectively. 
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Figure  2:  Photoelectron  spectrum  of  triazolyl:  upper  panel  experimental  spectrum  and  lower 
panels  computed  spectrum,  left  panel  minimum  energy  crossing  point  as  the  origin  of  Hd,  right 
panel  ground  state  minimum  as  origin  of  Hd, 

While  additional  work  on  the  representation  of  Hd  is  required,  the  improvement  when  the 
minimum  energy  conical  intersection  point  is  used  as  the  origin  is  palpable  and  reflects 
positively  on  the  weighted  normal  equations  approach  we  developed  in  the  course  of  our  study  of 
pyrrolyl  described  below. 

(iii)  The  photoelectron  spectrum  of  pyrrolide 

The  photoelectron  detachment  spectrum  of  pyrrolide,  the  anion  of  pyrrolyl  (CH)4N,  has 
been  the  object  of  recent  experimental  and  theoretical  studies. n  The  photodetachment  of  an 
electron  from  the  ground  state  of  pyrrolide,  produces  a  neutral  in  either  its  nominal  ground  2A2 
state  or  its  excited  2B,  state.  Here  the  attribute  nominal  is  meant  to  indicate  that  both  the  2A, 
state  and  the  “B,  state  are  actually  on  the  ground  state  potential  energy  surface  as  the  result  of 
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coupling  via  a  conical  intersection.  Interestingly,  the  experimental  photoelectron  spectrum  is 
devoid  of  features  attributable  to  the  :B,  state  and  those  features  resulting  from  the  nominal  :A: 
state  are  perturbed.  There  are  two  potential  sources  for  the  absence  of  spectral  features 
attributable  to  the  Bj  state,  strong  nonadiabatic  effects  or  small  transition  moments  (technically 
small  photodetachment  cross  sections).  This  diminution  of  the  intensity  of  peaks  in  an  adiabatic 
states  spectrum  by  nonadiabatic  and  transition  moment  effects  is  a  problem  of  enduring 
interest.14,12  This  issue  was  the  principal  concern  of  our  study. 

Attempts  to  describe  this  photoelectron  spectrum  using  an  adiabatic  state  model  were 
largely  unsuccessful.  The  only  existing  theoretical  nonadiabatic  simulation  of  this  spectrum 
used  a  simple  computational  model,  only  linear  as  opposed  to  fully  quadratic  vibronic  coupling 
was  included  and  only  a  limited  number  of  normal  modes  were  explicitly  treated.13  While  that 
analysis  did  demonstrate  the  importance  of  nonadiabatic  effects  in  the  diminution  of  the  2B, 
contribution  to  this  spectrum,  questions  concerning  the  origin  of  the  diminution  still  remained. 

We  have  completed  our  computational  study  of  the  photoelectron  spectrum  of  pyrrolide 
and  are  preparing  the  results  for  publication.  A  very  reliable  Hd  based  on  the  fully  quadratic 
vibronic  coupling  model  developed  using  weighted  data  averaging,  has  enabled  more 
quantitative  assessment  of  the  role  of  nonadiabatic  effects  in  this  photoelectron  spectrum  than 
was  possible  in  the  previous13  study.  This  Hd,  which  has  as  its  origin  at  the  minimum  energy 
point  of  conical  intersection  of  the  "A2  and  the  Bj  states,  describes  both  this  crossing  region  and 
the  minimum  energy  structures  of  these  two  states  very  accurately. 

The  simulated  and  experimental  photoelectron  spectra  are  presented  below.  The 
agreement  is  seen  to  be  excellent. 
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Figure  3:  Simulated  (  upper)  and  experimental2  (lower)  photoelectron  spectrum  of  pyrrolide. 


The  quality  of  the  simulated  spectrum  allows  us  to  reach  reliable  conclusions  concerning  the 
relative  importance  of  the  two  factors,  nonadiabatic  interactions  and  photodetachment  cross 
sections,  in  the  diminution  of  the  2B,  state  contribution  of  the  photoelectron  spectrum.  From  our 
study  it  emerged  that  the  intensities  of  the  simulated  spectrum  reflect  both  the  influence 
nonadiabatic  interactions  and  the  intrinsic  photodetachment  cross  sections  for  the  states  involved. 
There  are  currently  no  direct  measurements  or  calculations  of  these  photodetachment  cross 
sections  Flowever  our  careful  treatment  of  nonadiabatic  effects  allowed  us  to  draw  conclusions 
about  the  size  of  these  cross  sections  which  will  serve  as  invaluable  benchmarks  when,  in  the 
future,  calculations  of  the  photodetachment  cross  sections  become  available. 

D.  Photodissociation  of  nitramide  H2N-N02. 
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The  mechanism  the  photodissociation  of  high  energy  materials  such  as  HMX  and  RDX  is 
a  significant  issue  because  of  the  potential  role  of  excited  electronic  states  in  the  detonation 
process.8,9  Experimental  measurements  have  established8  that  the  initial  step  in  the 
photodissociation  of  these  and  related  nitramines  including  (CH3)2N-N02  is  the  production  of 
NO  following  a  nitro-nitrite  rearrangement.  Previous  computational  efforts  to  explain  NO 
production  did  not  focus  on  the  potential  role  of  conical  intersections  in  this  process  9  We 
believe,  that  based  on  other  research, this  neglect  is  not  justified. 

As  a  first  step  in  a  systematic  study  of  the  initial  step  of  the  photodissociation  of 
nitramines  vve  are  investigating  the  topography  of  the  S0  and  S,  potential  energy  surfaces  of 
nitramide.  The  results  of  our  first  order  configuration  interaction  treatment  based  on  a  DZP  basis 
are  summarized  in  Figure  4,  below. 

40000 


Figure  4.  Key  points  on  the  S0and  Sx  potential  energy  surfaces  of  nitramide  (filled  diamonds) 
connected  by  linear  synchronous  transit  paths  (open  diamonds)  on  S,  before  the  conical 
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intersection  (labeled  Intersection)  and  S0  after  conical  intersection.  The  abscissa  scale  is 
arbitrary.  Product  =  rearranged  structure  H2N-ONO.  Dissociated  product  =  H2NO  +  NO. 

This  figure  shows  that  following  excitation  to  S1  from  the  equilibrium  geometry  of  H:N- 
N02  on  S0,  the  molecule  can  evolve  through  the  region  of  a  minimum  on  Sj  and  over  a  barrier,  to 
an  Sq  —  Sj  conical  intersection.  To  consider  the  impact  of  passage  through  this  conical 
intersection  on  the  photodissociation  we  examined  “routing”,  the  motions  favored  by  the 
branching16  or  g-h  plane  as  the  wave  packet  emerges  on  the  ground  state  surface.  By  following 
the  gradient  from  distinct  points  in  the  g-h  plane,  we  found  that  after  going  through  the 
intersection  to  ground  state,  the  molecule  is  either  routed  back  to  the  ground  state  minimum,  or 
to  the  rearranged  product  H2N-ONO  which  subsequently  dissociates  into  NO  and  H2NO.  See 
Figure  4.  This  analysis  supports  our  conjecture  that  the  S0  -  S,  conical  intersection  seam  plays  a 
key  role  in  the  photochemical  production  of  NO. 

However  this  analysis  may  be  too  simplistic.  Although  only  a  single  point  of  conical 
intersection  (the  minimum  energy  crossing  point)  is  noted  on  Fig.  4,  an  extended  seam  of  conical 
intersection  near  the  nitro-nitrite  rearrangement  path  was  found.  This  continuous  seam  includes 
both  geometries  that  are  close  to  the  Sj  minimum  and  those  that  are  more  similar  to  the 
rearranged  product  H2N-ONO.  We  are  currently  analyzing  the  effect  of  routing  at  different 
points  on  the  S()  -  S,  seam.  In  the  future  we  will  extent  of  our  electronic  structure  treatment  to 
the  more  reliable  second  order  configuration  interaction  level.  Subsequently  more  complicated 
derivatives  of  H2N-N02  will  be  considered. 

E.  Future  Directions 

In  the  renewal  of  this  grant  we  intend  to  pursue  the  following  research  directions 
described  in  detail  in  our  renewal  proposal.  Determine  the  photoelectron  spectra  of  the  azolides, 
including  triazolide  ((CH)2N3'  -  work  in  progress)  and  pentazolide  (N5  work  just  stalled).  In 
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order  to  efficiently  perform  the  electronic  structure  calculations  required  to  determine  the  quasi- 
diabatic  Hamiltonian  Hd  for  penazolyl,  algorithms  previously  developed  to  efficiently  handle  the 
C3v  point  group  will  be  extended  to  arbitrary  nonabelian  point  groups.  We  anticipate  introducing 
Hd  with  higher  order  (  that  is  beyond  quadratic)  terms  and  ab  initio  calculated  transition  dipole 
moments  into  our  simulations  of  photoelectron  spectra.  Finally  we  will  consider  the  role  of 
conical  intersections  in  decomposition  of  electronically  excited  nitramines  and  other  energetic 
species  including  diazirine  H2CN2.  This  later  project  will  be  performed  in  collaboration  with 
Hanna  Reisler  at  the  University  of  Southern  California. 
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